Thermodynamics of 



Superfluidity 



A.F. Andreev and L.A 



Melnikovsky 



Abstract 



New, superfluid specific additive integral of motion is found. This 
facilitates investigation of general thermodynamic equilibrium conditions 
for superfluid. The analysis is performed in an extended space of ther- 
modynamic variables containing (along with the usual thermodynamic 
coordinates such as pressure and temperature) superfluid velocity and 
momentum density. The equilibrium stability conditions lead to thermo- 
dynamic inequalities which replace the Landau superfluidity criterion at 
finite temperatures. 
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Existence of two independent macroscopic motions is a specific feature of su- 
pcrfluids. Equilibrium state is, thus, characterized by two velocities of these 
motions. In general the equilibrium thermodynamic state is unambiguously de- 
termined by the values of all additive integrals of motion. Integral of motion 
that corresponds to one of the velocities (the "normal component velocity"), 
like it is in usual hydrodynamics, is momentum. 

Appearance of the superfluid velocity as a thermodynamic variable implies 
the existence of an extra additive integral of motion. Below we explicitly find 
the integral and this allows us to investigate general thermodynamic stability 
conditions and, in particular, the problem of critical velocities. 

Usually in experiments the vortices destroy superfluidity at Feynman critical 
velocity which is far below the Landau critical velocity. This is why the super- 
fluid hydrodynamics equations can be expanded in powers of low velocities and 
one safely uses the first nontrivial terms of this expansion. 

Nevertheless, there is a number of experiments 0E] where this is not enough. 
In some of them the superfluid flow is studied in small orifices. It has been shown 
that in these circumstances the maximum velocity is inversely proportional to 
the orifice width and may reach the order of the Landau critical velocity if 
the aperture is small enough. Other experiments are performed at tempera- 
tures extremely close to the superfluid transition. In this region, the linear 
approximation fails at low velocities starting from those much smaller than the 
Feynman critical velocity (which does not depend on the temperature). This 
means that all thermodynamic quantities of the superfluid become nontrivial 
functions of the not small superfluid velocity (i.e., they depend not only on the 
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usual thermodynamic coordinates such as pressure and temperature) . The only 
assumption one can make (and we do it) is that the fluid at rest is isotropic. 
This quite general statement of the problem is used in the paper (and in our 
earlier letter 3 ) ; in the next section we find the complete set of thermodynamic 
inequalities in this light, i.e., the conditions imposed on thermodynamic func- 
tions for the superfluid to remain stable. 

Later we employ the Landau phonon-roton model to calculate the highest 
velocity compatible with obtained thermodynamic inequalities and show that it 
can be interpreted as a critical velocity. This thermodynamic scenario suppos- 
edly explains the superfluidity breakdown in small orifices. In the vicinity of 
the A-point we use scaling speculations to find critical behavior of the critical 
velocity. 



Thermodynamic equilibrium of an isolated system is a state with maximum 
entropy at given values of all its additive integrals of motion. 

Superfluid system is unusual in the sense it has a continuous gauge symmetry 
spontaneously broken, such system possesses a specific quantum order param- 
eter, "wave function of the BE condensate". This is why in superfluid there 
are two independent velocities of normal and superfluid motion. The superfluid 
velocity is (up to a factor) a gradient of the wave function phase v s = V<f>. 
The same is true for the time derivative v s = V0. This is a local conservation 
law of all three components of the vector V s = J v s dr. In other words, v s is 
the density of a new additive integral of motion (along with mass density or 
momentum density). This integral is specific to superfluid. 

Instead of the entropy maximality, it is usually convenient 0] to use another 
requirement, equivalent to the first one, namely the energy minimality under 
constant entropy and additive integrals of motion. This is the requirement 
of thermodynamic stability and may be expressed in form of thermodynamic 
inequalities. 

Total energy of the superfluid E to t is an integral of the energy density E 
over the entire volume E to t = J E dr. The energy density can be obtained via 
a Galilean transformation[5] 



Here, v s is the superfluid velocity, p is the mass density and subscript "0" de- 
notes quantities measured in the frame of reference of the superfluid component 
(that is the frame, where the superfluid velocity is zero). Namely, Eq and jo are 
the energy density and the momentum density with respect to the superfluid 
component. The former is a function of p, jo, and the entropy density S. Its 
differential can be written as 
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E= ^+v s j +£ . 



(1) 



di?o = T dS + p dp + w djo , 



(2) 
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where Lagrange multipliers T, p, and w are the temperature, the chemical po- 
tential, and the so-called relative velocity of normal and supcrfluid components. 

The liquid is isotropic and, consequently, the velocity w and the momentum 
density jo are parallel to each other: 

w 

jo = jo(T,p,w)—. 

w 

This leads to a useful identity for the partial derivatives of jo with respect to 
w: 

8jq\ w k w l ( dj \ ( S kl w k w l 



dw 1 J T w 2 \ dw j T 



, >• (3) 
Further transforming QJ, we can write dE with the help of J2J) in the form 



c\E = T dS + I M + y - v s v n J dp + (j - pv n ) dv s + v n dj, (4) 

where j = pv s +jo is the total momentum density and v n = v s +w is the normal 
velocity. 

Stability implies that each "allowed" fluctuation increases the total energy 
of the system -Etot- Allowed are the fluctuations leaving conserved quantities 
unchanged. This means that the minimality of -Etot must be investigated un- 
der fixed entropy and all additive integrals of motion: mass, momentum, and 
superfluid velocity. 

Consider a macroscopic fluctuation of all the variables SS, Sp, <5v s , and 5j. 
They are conserved and this ensures that the first variation of the total energy 
for a uniform system is identically zero 

The matrix of the second variation quadratic form is a Jacobian matrix 8x8 
d(T, v n , p + vl/2- v s v n , j-pv n ) 



si 



9{S, j, p, v s ) 



It is positive definite if all principal minors in the top-left corner Mi, M2, . . . Ms 
are positive. We recursively test these minors: 

• The first positivity condition 

d(T, j, p, v„) _ d(T, j, p, Vs) d(T, v n , p, v s ) 



Mi 



d(S, j, p, v s ) d(T, v n , p, v s ) d{S, j, p, v g ) 

9 jo \ ( dj x 2 




dwJ Tp \dT 



> 



3 



corresponds to the usual requirement of the heat capacity positivity. It is 
shown below that (djo/dw) T > 0, hence the last inequality eventually 
becomes 



8S_ 
df 



p.W 



dw 



T,p 



djp 
dT 



> 0. 



(5) 



p , w 



• Positivity of the next group of minors is easily verified with the following 
transformation 



Si' 



d{T, v u , p, v s ) 
d(S, j, p, v s ) 



d(T, j, p, v s ) 
d{S, j, p, v s ) 



d{T, v n , p, v s ) 
d(T, j, p, v s ) 



(6) 



Positivity of the minors M2, M3, M4 is guaranteed if all minors of the 
second term in (JSJ are positive 



<9v n 



T,p,v s 



-1 


(8h\ 




\ dw /T,p 



(dj /dw) T<p 







jo/w 
jo/w 



Here we used (0 and chose the direction of the w vector as the first 
coordinate. Our collection is therefore augmented by two more inequalities 



jow > 0, 



djp 
dw 



> 0. 



(7) 
(8) 



T,p 



• The same transformation applied to the biggest minors yields: 



f>=& 



d{T, v n , fi + vl/2 - v s v n , j - pv n ) 



d(T, v n , 



p, 



) 



- si' Si' 1 



Again, the minors M5, Mq, M7, and Mg are positive if all nontrivial prin- 
cipal minors of Sj" are positive. We use the thermodynamic identity to 
relate the chemical potential p and the conventional pressure p 



dp 

This gives 

( d(n + v*/2-v B v n y 
[ dp 



dp 
P 



-dT- J ^dw, 
P P 



T,v n ,v s 



dp- 

9 PJt,; 
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The following is an explicit representation of $j" sub-matrix corresponding 
to a four-dimensional sub-space p, Vg,v^,Vg] as before we let the a;-axis run 
along w direction. Using © we obtain 



(dp/dp) T jp (djo/dp) q 



(dj /dp) Tw 





w 



(djo/dw) 








P - Jo/w 








P - Jo/w 



Appropriate inequalities are: 



T,w 



(9) 



which is literally a generalized (to a non-zero inter-component velocity w) 
positive compressibility requirement, 



Jo < wp, 



(10) 



dp\ / / dj \ \ (f dj 



Inequalities JSJ, Q, ©, ©, (fTOf) . and (fTTf) are sufficient conditions for the 
thermodynamic stability. 



3 DISCUSSION 

In a "stopped- normal-component" arrangement, the mass flux f with respect 
to the normal component may become more convenient than j — the mass flux 
relative to the superfluid one. The obvious relation between them / = pw — jo 
leads to the following reformulation of the inequalities: 

fw < 0, / < wp, (12) 

< < P- (13) 



5?U'-(£)J>(£L 

dP^T.w n,T > 9 i. d p\o, T ' ^ 
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As a simple application of the derived inequalities, consider them at w = 0. 
From O, JT3J), (JIU, and JTBJ) we get 

dS\ (d P \ 
„>(§*>) >o. 

Using conventional notation, last inequality reads 

Ps > 0, p n > 0. (16) 



4 PHONON-ROTON MODEL 

Here, we provide a practical example of the stability criteria for real superfluid 
4 He. To calculate derivatives involved in the inequalities one must use the 
microscopic approach. Simple and clear Landau phonon-roton model works 
pretty well over wide temperature and velocity ranges. We use this model to 
calculate the contribution of these quasiparticles to the "modified" free energy 
in the frame of reference of the superfluid component: 

Tq=Tq- wj , AT = ~S AT - j dw. 

This potential is obtained from the excitation spectrum with the conventional 
formula 

/ pw — e(p) \ \ dp 



T = T j In (l - exp 



T ) ) (27rft) 3 ' 

We denoted the excitation energy e(p), which is given for two branches by the 
expressions 

e P hb) = cp, e r (p) = A + (p - p Q ) 2 /{2m). 

Subscripts distinguish the quantities related to phonons and rotons, c is the 
sound velocity, A is the roton energy gap, m is the effective mass, and po is the 
momentum at the roton minimum (we take numerical data from Refs. jS], [7]: 
p = 0.145 g/cm 3 , A = 8.7 K, m = 0.16to H c, poP" 173 ^" 1 = 3.673 • 10 8 g" 1 / 3 , c = 
238 m/s, dA/dp = -0.47 • 10" 14 cm 5 S - 2 , dm/ dp = -0.45 • 10~ 23 cm 3 , dc/dp = 
467- 10 3 cm 4 s~ 1 g~ 1 ). A small dimensionless parameter ml\jp\ ~ 0.03 <C 1 
ensures, e.g., that the Landau critical velocity is determined by = A/pq. 

When integrated, these dispersion laws give the following contribution to 
the free energy: 

One can obtain all thermodynamic variables by differentiating this potential. 
Namely 

S = -(d? /dT) , j = -(d? /dw) . 

\ /w,p V / T,p 
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Figure 1: Critical velocity w c versus temperature T. Dashed line corresponds 
to the equation T = A — pqw. Note that the condition T < A — pow holds true 
over entire stability domain. The "stability" critical velocity w c coincides with 
the Landau critical velocity i>l at zero temperature and vanishes completely 
(see at the critical temperature T\. The curve is superimposed on the 

experimental graph of the critical phase difference over a 2 /im x 2 /jm aperture 
in a thin foil at various frequencies. [5j 



We neglect the quasiparticle contribution to the pressure derivative. 

At zero velocity the inequalities in Eq. fail at the A-point. When 

the velocity is higher, the inequality is the first to become invalid. At 
zero temperature, the critical (the highest possible without stability violation) 
velocity tends to the Landau critical velocity . Entire stability region is 
plotted in Fig. ^ an d Fig. [5] The liquid is unstable above the curves. 

In Fig. the critical velocity is plotted with the solid line. Scaled experi- 
mental datajS] are presented in the same graph. High frequency points follow 
the shape of the theoretical curve, but at low frequencies sufficient discrepancy 
is observed. This might be explained by the long time required for vortex nucle- 
ation. At higher frequencies the thermodynamic stability limit is observed while 
lower frequency critical velocity is cut off by vortices. Whatever feasible expla- 
nation this might be, the following warning must be taken into account when 
experimental data on critical phase difference and our predictions are compared. 
Assumption that the critical velocity is determined by the stability limit implies 
that the hydrodynamic equations inside orifice are essentially nonlinear. Par- 

^Note that for systems where all quasiparticles can be described hydrodynamically (in 
other words, systems without roton branch) inequality 1111 at zero temperature includes a 
condition (dp/dp) T — w 2 > 0, i.e., w < c. 
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ticularly, one can not consider the superfluid component as an incompressible 
fraction of the liquid. In other words the phase shift across the length of the 
orifice is not proportional to the maximum velocity attained in the orifice. 

In Fig.|21the critical energy flux is plotted. In typical arrangement the energy 
is carried by the superfluid counterflow, where normal and superfluid velocities 
are directed oppositely to each other to keep zero net mass flux. In zero mass 
flux frame the condition j = yields 

Q = Ts(w-^\ + (w- J -^j 30 (17) 

for the energy flux. 




T/T x 

Figure 2: Critical energy flux (at zero mass flux)<3 c versus temperature T. 



5 CRITICAL BEHAVIOR 

To find the behavior of the critical velocity near the A-point, scaling arguments 
should be used. We assume that the low velocity expansion of the irregular part 
of thermodynamic potential near T\ starts with 

Ar 2 - a + p s V l~Br x vi, (18) 

where A and B are constants, r = (T\ — T)/T\, a is the critical index of the 
specific heat C p oc r _Q , the "superfluid density" p s follows Josephson scaling 
relation 9 p s oc r^ 2 ""^ 3 , and the index x will be determined below. According 
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to the hypothesis of scale invariance, all terms in (|18fl are of the same order 
in the fluctuation region. This yields t x oc p 2 T Q ~ 2 and x = (a — 2)/3. On 
the other hand, the stability is destroyed by the velocity at which the linear 
approximation fails, i.e., 

W C OC y/ps/T" OC T (2 -" )/3 . (19) 

Current value for the specific heat index[TU| a = —0.0127 ± 0.0003 gives 0.671 
exponent for the critical velocity 

Extensive data near the A-point are available in the critical energy flux 
experiments. 2 In the vicinity of the transition, Eq. (|17|) is reduced to 
Q ~ TS(w — jo/ p) and the critical flux is determined by 

Qcocp^ocrt 4 - 2 ")/ 3 ^ 1 - 342 , (20) 

which is in moderate agreement with experimental value of the exponent 1.23 ± 
0.02. 

6 CONCLUSION 

Experimentally, the superfluidity breakdown is believed to have the following 
nature. J] Until the aperture size is too small, the critical velocity does not 
depend on the temperature and increases as the size decreases. This is the very 
behavior that is specific to the Feynman vortex-related critical velocity. 

When the orifice width is narrow enough, or the temperature is close enough 
to T\ the breakdown scenario and its features change. The critical velocity does 
not depend on the aperture any more but decreases when the temperature 
increases. This behavior is commonly associated^ with the Iordanski-Langer- 
Fisher mechanism. ( 1 1 j Nevertheless, this association lacks numerical comparison 
because no reliable information about the actual orifice shape is available. 

On the other hand experimentally observed behavior of the critical velocity 
can be attributed to the suggested stability criterion. In other words, we provide 
an alternative explanation of experimental results based on an assumption that 
in narrow orifices the thermodynamic limit of w c is reached. 

We should also note that our approach to the critical velocity as a stability 
limit is similar to that used by Kramer. ^21 Actually, the inequality he employed 
is not a thermodynamic one, but numerical results for the critical velocity he 
obtained using the phonon-roton model do not deviate much from those plotted 
in Fig. [2 
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